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Abstract 

Ecologists and conservation biologists need to identify the relative importance of species to make sound 
management decisions and effectively allocate scarce resources. We introduce a new method, termed env- 
iron centrality, to determine the relative importance of a species in an ecosystem network with respect to 
ecosystem energy-matter exchange. We demonstrate the uniqueness of environ centrality by comparing it 
to other common centrality metrics and then show its ecological significance. Specifically, we tested two 
hypotheses on a set of 50 empirically-based ecosystem network models. The first concerned the distribution 
of centrality in the community. We hypothesized that the functional importance of species would tend to be 
concentrated into a few dominant species followed by a group of species with lower, more even importance 
as is often seen in dominance-diversity curves. Second, we tested the systems ecology hypothesis that indi- 
rect relationships homogenize the functional importance of species in ecosystems. Our results support both 
hypotheses and highlight the importance of detritus and nutrient recyclers such as fungi and bacteria in 
generating the energy-matter flow in ecosystems. Our homogenization results suggest that indirect effects 
are important in part because they tend to even the importance of species in ecosystems. A core contribution 
of this work is that it creates a formal, mathematical method to quantify the importance species play in 
generating ecosystem activity by integrating direct, indirect, and boundary effects in ecological systems. 

Key words: ecological network analysis, network environ analysis, food web, trophic dynamics, Ecopath, 
functional diversity, biodiversity, ecosystem function, foundational species 



1. Introduction 



Identifying the functional importance of species is a fundamental step in describing community and 
ecosystem structure and function. It is essential for ecologists, conservation biologists, and resource man- 
agers to understand species relative importance so that they can make infor med management decisions and 
effectively allocate limited conservation resources. Rank-abundance curves (jWhittakerl . I1965T ) are a classic 
tool used to describe the biodiversity and relative importance of species in a community. This approach gen- 
erally assumes that the importance of the species is proportional to its abundance or biomass. Alternative 
techn iques include using functional measures like productivity i n dominance-diversity c urves (jWhittakerl . 
Il965h. and more rec ent concepts such as ne twork ro l e equivalence dLuczkovich et all l2003h , keystone species 



( Power et al . 19961) . ecosys tem engineers ( Lawtonl 1994 J° nes et al, . 1994 ). and more generally the con- 



cept of foundational species (|Davtonl . [T973lEllison et al.l . l2005l) . However, quantifying the relative functional 
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importance of species embedded in their communities remains a challenging problem. In this paper, we in- 
troduce a new approach to quantifying the relative functional importance of species in ecosystem networks 
termed environ centrality. While species importance can be defined in a variety of ways, here we specifically 
focus on their importance for energy-matter distribution in communities. While we focus on ecosystems, 
the methods are generally applicable to any complex system of conservative fluxes modeled as a network. 

One reason that the fundamental step of quantifying species functional importance has not been fully 
addressed is that indirect interactions complicate the assessment. Organisms are embedded in an intri- 
cate web of interactions and it is this tangled network of relationships that en ables indirect influences 



*,o become signilicuni couiponnits of ecolo»i<-i>l and ovolulioiutry interact i< his fe.g. Patten. IttsJ: Wool Ion 



19941: IWhitham et all 1 20061 iBascompte and Jordanol l2007t lEstrada and Bodinl . 120081: iMontova et all 12009; 



Keith et alll2010h . For example. IPoulin et al.l(|2010l ) showed how invasive species can increase the prevalence 



and severity of disease in a community through trait-mediated indirect effects. Troph ic cascades are one type 



of trophic indirect ef fect that can h ave large and unexp ected impacts on ecosystems (|Carpenter et all 11985 



Moonev et al. . 2010f) . For example, Berger et al 



2008T) found that where wolf populations were extripated, 
mesopredator populations like coyotes were released from predator control. This then resulted in a four fold 
decrease in the survival rate of pronghorn fawns, a coyote prey item. Thus, understanding the ecological and 
functional importance of species in their ecosystem requires understanding the full environment of direct 
and indirect influences. 

Network models of ecosystems let ecologists quantify species importance and are a key tool for deter- 
mining the importance of indirect influences that are distributed through these ty pes of interactions (e.g. 
iBondavalli and Ulanowicz , Il999t iDame and Christianl l2006t Bclg rano et al.1 1 120051) . In these models, net- 
work nodes represent species, groups of species, or abiotic resource pools. Although the nodes may be 
a species complex or non-living compartment of energy-matter, we generically refer to them as species 
throughout this paper for simplicity. Network links represent the flow of energy-matter from one node to 
another. These energy-matter flows can be created by several ecological processes including feeding, death, 
and excretion. This representation of complex ecosystem interactions lets ecologists apply mathematical 
and computational tools to learn more about the structure and function of ecosystems. These ecosystem 
networks can be represented by weighted and directed graphs so that a link not only implies a relation- 
ship between two nodes but it also indicates how much (weighted) energy-matter and from whom to whom 
(directed) the energy-matter is flowing. 

Ecological Network Analyses (ENA) are a collection of quantitative methods that systematically asses 
information from a full, co mplex network description. Th ere are several specific imple mentations of this 
concept, such as Ecopat h ( Christensen and Walters . 2004) . Ne twork Environ A nalysis (Fath and Patten . 



concept, such as tcopath (Christensen and Walters, iNetwork Environ Analysis (.bath and rattei 

1999bl : iFath and Borrettl l2006h . EcoNet (jSchramski et al.L l201oh . and NETWRK (|Ulanowiczlll986l) . ENA" 



throughflow analysis mathemat ically partitions e c osystem energy-matter flow acr o ss boundary, direct, an d 



indirect paths in the network ( Fath and Patter] . Il999bl : Fath and Borrett . l2006t ISchramski etail l2010h 



This a pproach has been applied to analyze the structure and function of ecosystems. For example, Finn! 
1 1980 ) used the tech niqu e to characterize minera l cycling in the Hubbard Brook ecosystem. More recently, 



Borrett et al.l (120061) an d ISchramski et al.l ( 2006 ) investigated nitrogen cycling in the Neuse River Estuary, 



and Zhang et al. ( 2010h investigated the urban water metabolic system of Beijing. ENA lets scientist study 
the effere nt and afferen t holistic environments of species within a system boundary, which Patten termed 



environs (|Pattenl . 119781 ). 



Based on the development and application of ENA, system s ecologists have hypothesized that indirect 
energ y -matter flow s tend to dominate direct flows in ecosystems (jHigashi and Patten , 1989; iJprgensen et all 
120071 ) . iFathl (|2004l) found support for th i s hypo thesis in large ecosystem models generated from a community 
assembly type algorithm, iBorrett et al.l (|2006l ) found supp orting evidence for the hy pothesis in 16 seasonal 
nitrogen cycling models for the Neuse River estuary, and Salas and Borrett] ( 2011 ) found general support 
for the hypothesis in 50 empirically-based trophic ecosystem models. A consequence of t hese large indirect 
energy-matter flows is that reso urces tend to be more evenly distributed in the systems ([Fath and Patten . 

>99a[ IBorrett and Salasi l2010h . Given this tendency, we hypothesized that indirect effects tend to ho- 
mogenize the relative importance of the species, decreasing the relative influence any single species has on 
ecosystem functioning. 
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This paper has three main objectives. First, we introduce a new metric of functional importance based on 
the throughfiow analysis of ENA and centrality concepts from social science. We contrast this measure with 
other existing centrality indexes to demonstrate its utility and uniqueness. Second, we use this measure 
to characterize the relative importance of species in 50 trophic ecosystem models. Third, we test two 
hypotheses regarding the functional organization of ecosystems. Based on previous community analysis, we 
first hypothesized that there will tend to be a "concentration of dominants" or functionally more important 
species and a "long tail" of species with a lower and more similar importance. We expected that the 
ecosystem dominants would be compartments like detritus, particulate organic matter, and bacte ria because 
of their critical role i n reci rculating energy-matter and connecting the green and brown food webs (jWilkinsonl . 
20061 : Jordan et al. . 2007 ). We further tested the hypothesis that indirect effects tend to homogenize the 



relative importance of species. 

This paper is organized as follows. In the next section we introduce the centrality concept and describe 
its use and diversity. We then describe the ENA throughfiow analysis that we use and introduce the environ 
centrality measure. Section 4 describes our analytical approach to testing our hypotheses and Section 5 
describes our results. We conclude the paper by putting these results into a broader context. 



2. Centrality 



Centrality is a concept used by scientists studying complex systems that was initially developed by so 



92 rial scientists t o quantify the importance of individuals in network models (Wasse rman and Faustl . [1994; 
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Borgattil 120051 ) . Metrics based on this concept indicate how a node's position in the network contributes 
to its st ructural or functional i mportance. There are mul t iple c entrality metrics with varying charac- 



teristics^ Wasserman and Faust, 1994; Borgatti and Everettl . [2006), but many tend to be well correlated 
(jValente et all l2008t iNewmanl . I2006TI . Despite this correlation, differences between the measures can be 
useful. For example, degree centrality is based on the number of network edges directly touching a node 
and describes the local or proximate importance of the node. In contrast, eigenvector centrality describes 
the stable distribution of pathways (i.e. at long pathways) passing through the nodes, which provides a 
more global or whole network understanding of the nodes' relative importance. These different measures 
have been well described in the literature, so here we focus on the ecological relevance of selected centrality 
measures. 



There are several approaches to quantitatively describing ecosystem networks (jUlanowiczi . 1986; Bc rsier et al 



2002) . Several metrics commonly used to desc ribe food webs such as link density and connectance (jDunne et al 



2002 ) are what IWasserman and Faust] (1994) call group level indicators of centrality. Link density is the 



average node degree in the network. Connectance is link density proportional to the size of the network. 
Degree centrality, and thus link density and connectance, is a local centrality metric because it only consid- 
ers the immediately adjacent neighborhood of connected nodes. Local centrality metrics may be a useful 
starting point for describing a food web, but they may not be the best descriptors of the organism's impor- 
tance in energy-matter transference because they assume the influence is restricted to a local neighborhood. 
Thus, they neglect the important indirect influences. From a trophic point of view, indirect effects become 
important in part because energy-matter is passed from one organism to another and may ult imately reenter 



the same organism throug h nutrient cycling, creating a well connected web of interactions (jAllesina et al 
20051 : iBorrett et all [2007n . 



Centrality metrics that incorporate indirect influences include betweenness centrality dFreeman , 1979 ) , 
eigenvector centrality ([Bonacich , ll972t lBorgattll2005tlEstrada and Bodinl. 120081 [ Allesina and Pascuali 120091), 
and weighted topological importance (jBauer et all l2010t IJordan et all 120031 ). These centrality metrics are 
more appropriate for many ecological applications than those that only consider direct influences, but still 
have important limitations. Betweenness centrality considers the possible gatekeeper role that some nodes 
may play in the transfer of information, energy, and matter from two more distinct groups. While th is metric 
can be useful, its focused on a particular type of bridging importance ([Wasserman and Faustl 1 1994) and was 
not intended to provide a more general metric of the importance of a node. In trophic or biogeochemical 
network models, eigenvector centrality effectively quantifies the equilibrium distribution of energy-matter 
flowing through the nodes, thus considering the direct and indirect interactions. However, there are two 
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potential problems with this measure for ecological networks. First, it ignores the contribution of the initial 
transient effects that maybe important in some contexts, especially systems in which the strength of interac- 



tion decreases rapidly with path length like ecosystems ( Borrett et al . 201(1) . Second, it only considers the 



dominant eigenvector, which dep ending on the structu r e of t he network may not be an adequate approxi- 
mation of the transfer dynamics ( Borgatti and Everett . 20061) . Weighted topological importance quantifies 
the effect a species has on others in a community, which is great for understanding competition but does 
not provide information on how species contribute to global network properties such as total energy-matter 
throughflow. To address these limitations, we introduce environ centrality, which uses weighted informa- 
tion, integrates direct and indirect effects, describes how species contribute to global network measures, and 
captures the transient dynamics as well as the equilibrium effects. 
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3. Ecological Network Analysis 

3.1. Throughflow Decomposition 

As comprehensive summaries of ENA methodology exist ( Ulanowiczl 1986 ; Fath and Borrett . 2006 



Schramski et al. . 2010f) . here we focus on the components of ENA necessary to calculate the central 



139 i ty of species in ecosyst em network models. We used the Network Environ Analysis Matlab® function 



i4o (jFath and Borrettl . 120061 ) to implement these analyses. 

hi Let a network model with n nodes be represented by a matrix, F nx „ = (fij), which defines the quantity 

142 of energy-matter being transferred from node j to node i. The structural component of F is the adjacency 

143 matrix, A = (ay), in which a 1 indicates a direct connection from j to i and indicates none. Energy- 

144 matter entering the system at node i is denoted by Zi , whereas energy-matter leaving the system at node j 
is y 3 . 

ENA uses this system data to determine several ecosystem properties. First, the total flow through a 
node is defined as throughflow, which can be input Ti = 2i+^j=i fij or output Tj — J27=i fij + Vj oriented. 
ENA typically assumes that the networks are at steady state (e.g., mass-balanced) so Xi = Tj. Second, we 
calculate the direct flow intensity matrix G. We focus on the output oriented direct flow intensity matrix, 
which is calculated as G = (gij) = fij/Tj. Elements of G are unitless and represent the direct flow intensity 
from node j to node i. Third, raising G to a power m, gives the flow intensity between two compartments 

152 over all paths of length m. The sum of flow intensities over all pathways in a network is defined as the 

153 integral flow matrix N, which is 



N=^G m = G° +GL + G 2 + ... + G t " + ... . (1) 

m — Boundary Direct Indirect 

154 N quantifies the intensity of output oriented throughflow from j to i over all pathways in the network. Since 

155 ecological systems are thcrmodynamically open, the energy-matter must dissipate as path length increases 

156 causing the sum of flow intensities over all pathways to converge. Thus, we can use the identity (I — G) _1 , 

157 where I = G° is the identity matrix, to find the exact values of N. 

158 Multiplying the integral matrix N by the input vector z returns the output oriented node throughflow 

159 T = Nz. This equation ensures that equation (TTJ) is a true partition of flow and that the flow elements 
wo are not double counted. Total System Throughflow (TST — Y^j=i ^i) ls the sum of node throughflows and 

161 is a global measure of the total network activity or size of the system. Thus, the integral matrix shows 

162 how TST is generated by species in an ecological network and incorporates energy-matter flux over all 

163 indirect pathways. Environ centrality is calculated from the integral flow matrix and quantifies the relative 

164 importance of species in creating total system activity. 

165 3.2. Environ Centrality 

166 Environ centrality is a measure based on the ENA output-oriented integral flow matrix, N. We introduce 

167 three related environ centrality measures: input (EC m ), output (EC out ), and an average of the two (AEC). 
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168 These are calculated as follows: 
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(2) 



EC m and EC out indicate the relative importance of species in generating ecosystem activity from an input 
or output direction, respectively. AEC determines the relative importance of species in ecosystem models by 
averaging the input and output importance. Thus, AEC quantifies the relative contribution of a node to the 
energy-matter exchange within a system. It is a global centrality measure of the functional importance of 
species because it incorporates direct and indirect pathways, weighting pathways by the amount of energy- 
matter flux intensity passing through a node. Further, it captures the transient dynamics as well as the 
equilibrium dynamics and describes how a species contributes to total s ystem activ i ty. F ig. [T] shows an 
example rank-order distribution of AEC for the Georges Bank ecosystem ( Link et all . 12008 ). 

ENA is an Input-Ou tput analysis in which it is possible to calculate output oriented and input oriented 



integral flow intensities (|Fath and Borrettl . 120061 : iSchramski et all l2010h . Thus, it is possible to calculate 
the EC m , EC out , and AEC for both of these orientations. The interpretation of output and input oriented 
integral flow intensities differ. The output orientation looks forward in time and follows the direction of the 
network arrows. The input orientation looks backward through the network, against the arrow directions, 
to see where the fl ow has come from. Both orientations can be analytica lly useful independently and when 
jointly considered ()Ulanowicz and Puccia , 19901 : Schramski et al. . 2006 ). As the purpose of this paper is 
to introduce environ centrality and illustrate its ecological significance, only the output oriented case is 
analyzed for simplicity. We expect the input oriented results to be similar, although the interpretation will 
differ. 



187 4. Experimental Design 

188 In this section we describe our experimental design to test our hypotheses. We first characterize the set 
is? of 50 empirically based ecosystem models we analyzed. We then explain our evaluation of the utility and 

190 uniqueness of EC. Finally, we describe our approach for determining species dominance and the centrality 

191 homogenization in 50 ecosystem networks. 

192 J^.l. Network Models 

193 We calculated AEC for 50 empirically-based network models of ecosystem energy-matter flows (Ta- 

194 ble[5]). iBorrett and Salasl ( 2010l) first assembled this data set to test the distinct systems ecology resource 

195 homogenization hypothesis. The collection is available from |http : / /people . uncw.edu /borretts /research. html| 
we Models were selected that spanned a range of sizes from smaller, highly aggregated systems comprising 4 

197 nodes to larger, less aggregated systems comprising 125 nodes. Each network model traces trophic and 

198 non-trophic energy-matter transactions such as feeding, excretion, mortality, immigration and emigration 

199 and together represent 35 distinct ecosystems. To avoid a selection bias, all models discovered in the lit- 

200 erature that wer e created to model a spe cific ecosystem and that included empirical data were included in 

201 the data set (see lSalas and Borrettl . 120 111 for more details). To meet the ENA steady- s tate a ssumption, we 

202 balanced 15 of the 50 models using the AVG2 algorithm, which Allesina and Bondavallil ( 2003 ) demonstrated 

203 minimized balancing bias. 

204 Jf.,2. AEC Sensitivity and Uniqueness 

205 Given the environ centrality metric, our first step was to evaluate its sensitivity and demonstrate its 

206 novelty. To illustrate its sensitivity, we compared AEC values for four realizations of a hypothetical ecosys- 

207 tern flow network (Fig. [2]). The four realizations vary in their distribution of flow amongst the five internal 
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direct pathways. This is evident in the direct flow intensity matrices, G, and the subsequent integral flow 
matrices, N. The direct flow intensity values in the second realization are 10% of those in the first, which is 
a universal extensive change. The distribution of flow intensities is the same, but in the second realization 
the total magnitude is less. Realizations three and four have equal total flow intensities, but differ in how 
it is distributed. Thus, this is an example of an intensive change. Together, these realizations demonstrate 
AECs ability to capture both intensive and extensive network changes. In addition to AEC , we also ex- 
amined the average eigenvector centrality AEVC (Table [T]) since these were likely to be the most similar 
ccntrality measures. 

To show AEC^s uniqueness, we compared its results in two ecosystem models to five commonly used cen- 
trality metrics: degree centrality (DC), weighted degree centrality (WD), weighted topological importance 
(WI), and betweenness centrality (BC). The Oyster Reef and Chesapeake Bay models were selected as case 
studies to demonstrate the distributions of these five centrality metrics in empirical network models (Fig. [3]). 
We calculated all metrics usin g common formulations (Table [I} and then normalized each set of scores by 
its sum (sensu Ruhnau . 2000fl . creating proportions to facilitate comparisons among the five metrics. In 



addition to the comparisons in Fig. |3l we compared AEC and AEVC for four network models; the Oyster 
Reef, Chesapeak Bay, Bothnian Bay, and Northern Benguela. Finally, we calculated the Spearman rank 
correlation of the AEC and AEVC in the 50 empirically based ecosystem models listed in Table O We 
focus on these two centrality measures because they are most likely to be similar because they are both 
global pathway centralities. 

4-3. Dominance and Evenness 

To address our hypothesis regarding the relative importance of species in ecosystems, we used the co- 
efficient of variation (CV) to characterize the evenness of the AEC distribution and to identify dominant 
species in 50 ecosystem models. Although there are multiple methods for quantifying variation, we chose CV 
because it is a non-dimensional measure normalized to the mean values. This lets us compare the relative 
variation across s ystems, even if they differ in t he number of species or TST . In addition, CV is sensitive to 



distribution skew (jFraterrigo and Rusakl . 120081 ) . which we use to our advantage to identify dominant species. 



We considered ecosystems with a CV less than unity (1) to be low variance. To identify the dominant 
species, we ranked species according to their AEC scores and calculated the CV for the entire community. 
We then progressively removed the highest ranking species until the CV of the remaining community AEC 
was less than or equal to 0.5. A CV value of 0.5 is low enough to identify dominant species without being 
so low as to claim all species are dominant. It also identifies which species are responsible for the highest 
proportion of the variation in EC. All species removed before reaching the CV threshold are classified as 
dominant species with respect to the ecosystem flux of energy-matter. Thus, this approach determines both 
the evenness of non-dominant species, and identifies the dominant species. 

4-4- Homogenization and Indirect Effects 

Our final hypothesis was that indirect flows would homogenize the relative importance of species in 
generating TST. We addressed this hypothesis by comparing the relative importance of species in the direct 
flows to those of the integral flows. Recall that equation Q] calculates the integral flow intensity matrix N 
that combines flow intensity over all boundary (G°), direct (G 1 ), and indirect pathways (J2m=2 G m ). To 
find the relative importance of species from the perspective of the direct flows alone, we substituted G for 
N in equation [2] to calculate direct average environ centrality (AEC direct ) . We then compared AEC dtrect 
to the integral average environ centrality (AEC) to test the homogenizing effect of indirect relationships on 
AEC. 

As the integral matrix includes boundary, direct, and indirect flow intensities, it is possible that observed 
homogenization could be caused by the boundary input as well as the indirect flows. To isolate the effect of 
the indirect flows, we also compared the AEC dlrect to the AEC calculated on N — G° instead of N. 

Homogenization of species importance was quantified by comparing the CV of AEC dlrect and the CV of 
AEC. We created a ratio of the two, CV(AEC direct )/CV(AEC), such that when the ratio is greater than 
unity (1) it indicates that the AEC is more evenly distributed than AEC dlrect . Ratio values greater than 
unity indicate that indirect effects homogenize the importance of species in generating ecosystem activity. 
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258 5. Results 



259 5.1. EC Sensitivity and Uniqueness 

260 To establish the sensitivity of AEC to both intensive and extensive changes, we applied it to four realiza- 

261 tions of the hypothetical ecosystem model shown in Fig. [5] The AEC distributions for the realizations are 

262 clearly different (Fig. HP). In the first realization, the detritus box was more important and the importance 

263 of the primary producer was diminished. In the second realization in which the direct flow intensities were 

264 reduced by 90%, the AEC values are much more similar. This is because the drop in flow intensities is 

265 then transmitted through the longer pathways, effectively discounting them. However, there is no difference 

266 between the AEVC centrality distributions between realization one and two, demonstrating that eigenvector 

267 centrality is not sensitive to this extensive change. 

268 Realizations three and four maintained the total magnitude of the network but have different distri- 

269 butions, which is an intensive change. The AEC and AEVC distributions between realization three and 

270 four are different, demonstrating that both metrics are sensitive to these intensive changes. In the third 

271 realization, one flux to the consumer compartment increased which caused it to become the second most 

272 central compartment. The detritus compartment also responded to the change, increasing its centrality 

273 from 0.26 to 0.30. For the fourth realization, the flux from detritivore to detritus increased, causing detritus 

274 to maintain its centrality from the previous realization at 0.30 and detritivores to increase its centrality to 

275 0.30 as well. These four cases show that AEC is sensitive enough to capture important differences between 

276 the model flow realizations, and specifically demonstrates its ability to capture both intensive and extensive 

277 changes to direct energy-matter flows. 

278 The comparison of the four centrality metrics to AEC in the Oyster reef and Chesapeake Bay ecosystem 

279 models (Fig. shows the different structural and functional importance of the species within each network 

280 and illustrates the unique perspective of AEC. DC highlights which species were most connected within the 

281 network topology (unweighted), yet was not always a predictor of which species would be most important 

282 in generating TST. There were several disparities between DC and AEC which show that being the most 

283 connected does not necessarily result in being a key species in generating TST . For example, the suspended 

284 particulate nutrient node from the Chesapeake Bay model has a greatly reduced ranking in AEC despite 

285 having the second largest DC. WD had similar mismatches with AEC, and demonstrates that high or low 

286 magnitudes of direct flows are not always indicative of overall importance in generating ecosystem activity. 

287 In both models, WI consistently placed primary producers and herbivores (i.e. filter feeders, phytoplankton, 

288 and zooplankton) as species with the highest effect on others, yet these species ranked low in AEC. The 

289 BC distribution was the most distinct centrality metric. 

290 Average eigenvector centrality should be the most similar to average environ centrality because they are 

291 both global centrality measures related to pathways or energy-matter flow; however, we see from Fig. [5] 

292 that they can differ. Fig. [4] further illustrates differences between AEVC and AEC. It shows both types 

293 of centralities for four cases: the Oyster Reef and Chesapeake Bay models used in Fig. [3] and the Bothnian 

294 Bay and Northern Benguela ecosystem models. The centralities indicate a difference of the rank-order 

295 importance of the nodes for all models except the Oyster Reef. For example, AEC suggests that Bacteria, 

296 DOM, Mesozooplankton, and Microzooplankton have a larger functional importance in the Bothnian Bay 

297 ecosystem than is suggested by the eigenvector centrality (Fig. Hp). In contrast, eigenvector centrality 

298 indicates a very dominant importance of Meiofauna, Sediment carbon, and Macrofauna. The two centralities 

299 also show a different rank importance of the nodes in the Northern Benguela ecosystem model. From the 

300 AEC perspective, POC has a dominant importance, while most of the other species have very similar 

301 centralities. Variation in AEVC is larger in all models, including Oyster Reef where the rank importance are 

302 the same. One interesting difference is the diminished importance of Seals suggested by AEVC for Fig. Hp). 

303 Additionally, AEVC ranks bacteria in sediment POC as more important than Sediment Particulate Carbon 

304 in Fig. UP while the AEC values are nearly the same. 

305 To generalize our analysis, we calculated the Spearman rank correlation between AEC and AEVC for 

306 the 50 ecosystem network models (Tabled]). These correlations ranged from 0.2 to 1 and had a median 
so? value of 0.92. As we expected, this suggests that these two centrality measures are typically, but not always, 
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308 similar. Our case studies in Fig. @] highlight how these small correlation differences can be ecologically 

309 meaningful. 

310 5.2. Dominance and Evenness 

3n Rank AEC curves provide a tool to visualize the relative importance of species (Fig. [TJ). The AEC 

312 distribution of the Georges Bank illustrates the tendency for a few nodes to have high AEC with most 

313 having relatively lower and even AEC values. AEC variation in the ecosystem models we analyzed was 

314 generally low (Table [3]). Twelve of the models have a CV(AEC) less than 0.5 and 30 of the models have 

315 a CV(AEC) between 0.5 and unity. However, eight of the models exhibit more variability with CV values 

316 larger than unity. In the models with a CV larger than 0.5, no more than three dominant species had to be 

317 removed before the CV of the remaining species fell below 0.5 (Table [3]) ■ The median number of dominants 

318 was two. As expected, detritus and detrital recyclers are predominantly responsible for generating TST in 

319 most ecosystems, with water flagellates being the only non-detrital or bacterial species to rank as dominant 

320 AEC contributors. 

321 The second hypothesis was that sub-dominant species would have a more even distribution of importance. 

322 In Fig. [TJ the first two species were considered dominant species and the CV score for the whole ecosystem 

323 was 0.68. However, once the dominants were removed the CV score of the remaining species was 0.41 

324 suggesting less variation in the importance of the remaining species. Table [3] identifies the CV value of 

325 the entire community, as well as the community without the dominant species for all 50 ecosystem models. 

326 After the dominant species were removed, the AEC scores were relatively even, with an average CV score of 

327 0.38. This indicates that the standard deviation was only 38% the magnitude of the mean. All the networks 

328 support our hypothesis that the importance of sub-dominant species would be relatively even. 

329 5.3. Homogenization and Indirect Effects 

330 The comparison of direct and integral AEC in the Oyster Reef model (Fig. [5]) illustrates how the relative 

331 importance of species in generating TST can become more similar when indirect flows are considered. Fig. [5] 

332 shows that this is a general trend as 49 of the 50 network models (98%) exhibited less variation in AEC 

333 when indirect flows were considered. The EMS estuary was the only network to not meet our expectation. 

334 Furthermore, when boundary effects were excluded from our analysis, 33 of 50 ecosystems (66%) exhibited 

335 homogenized importance due to indirect effects alone. Thus, indirect effects appear to homogenize the 

336 importance of species in ecosystems. 



337 6. Discussion 

338 6.1. Ecological Significance 

339 Our results make two ecologically significant contributions. First, AEC is a novel way to quantify the 

340 relative functional importance of species or functional groups in an ecosystem with respect to the global 

341 energy-matter exchange of a network. It is sensitive to both extensive and intensive changes, and incorpo- 

342 rates indirect effects. With this metric, we gain new insight into the distribution of species contribution to 

343 the total system activity. Second, the results highlight another critical consequence of indirect interactions 

344 in ecosystems; they tend to make species' functional importance more similar. If this pattern holds generally, 

345 it has important implications for conservation biology which we discuss below. 

346 Our results support the hypothesis that ecosystems are comprised of a few functionally dominant species 

347 and a larger set of species with relatively similar functional importance. The smaller CV values for the 

348 remaining community after the dominants are excluded demonstrate that the top ranked species are respon- 

349 sible for a majority of the variation in AEC and that most species in ecosystems have relatively low and 

350 even centrali ties. This is furth er evidence for the general dominance-diversity community patterns initially 

351 described bv lWhittaken (Il965h . 



352 The second hypothesis this paper tested concerned the effect of indirect flows on the species' environ 

353 centralities. Our results support the hypothesis that indirect effects tend to even the functional importance 

354 of species. Indirect flows alone homogenized the importance of species in 66% of the ecosystem networks 



355 analy zed. The empirical literature is rich with examples of indirect relationships driving ecosyst em reac- 

356 tions ( Berger et all 2008 ; Bondavalli and Ulanowicz , 1999t Carpenter et al. . 1985 ; Wootton , 1994 ) and our 

357 theoretical analysis provides evidence for another important consequence of indirect effects. 

358 The EMS estuary was the only network in which the integral centrality AEC was not more even than 

359 the direct environ centrality AEC Dlrect . This could be because of the extreme dominant centrality of the 

360 sediment POC in this model. After removing sediment POC from the AEC variation analysis, the CV 

361 declined 70% which was one of the largest percent drops in variation. It was also one of the few ecosystems 

362 to have a CV(AEC) that was greater than unity. 

363 Collectively, our homogenization results support the trend towards whole ecosystem management for 

364 conservation biology. Although a few species may have a dominant functional importance, most species are 

365 contributing in similar ways. This could be interpreted as suggesting that losing the non-dominant species 

366 would have less significant impacts on the system, but this misses a fundamental aspect of the results. The 

367 homogenization of the species centralities is a consequence of the indirect interactions that are distributed 

368 across the network of relationships. This depends on the existing species and their pattern of direct and 

369 indirect connections. Thus these results show that preserving the intact network is critical to maintain 

370 current ecosystem total energy-matter throughflow. 
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377 
378 
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380 
381 
382 
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384 
385 



389 
390 



6.2. Centrality metric comparisons 

Our comparison of centrality metrics in Fig. [3] highlights the differences and similarities between these 
five metrics. In the case of WD and WI, both metrics are calculated from the same matrix F, but differ by 
the anchoring metric (i.e. sum of energy-matter entering a node versus total energy-matter flowing in the 
network) and the length of the pathways considered. This may explain the similarities between WD and 
WI as seen in the Oyster Reef model. 

Average eigenvector centrality should be the most similar to average environ centrality because they are 
both global centrality measures related to pathways or flow. However, our results show that these centrality 
measures can be quite different. The results of the hypothetical ecosystem model are most instructive 
(Fig. The first and second model realizations whose direct flow intensity distributions are the same have 
the same AEVC distributions. The large drop in direct transfer efficiencies does not affect the equilibrium 
distribution of flows, thus the eigenvector centralities remain the same. This is unlike AEC which does 
respond to the drop in transfer efficiency because it captures the transient flow dynamics as well as the 
equilibrium distribution. This demonstrates the ability of AEC to capture extensive changes that AEVC 
does not. Both AEC and AEVC capture differences between realizations three and four, demonstrating 
that both are able to capture intensive changes to direct energy-mat ter flow. Although AEC and AEVC 



are well correlated, most centrality metrics tend to be well correlated (jValente et all l2008t iNewmanl . [2006) 



and small differences between metrics can be biologically important (see Fig. @|. 

Alternative centrality measures have different values. We do not expect one to be the universally correct 
analytical choice. If the goal is to understand the relative importance of nodes in generating the total system 
activity of a network model of conservative flows, the evidence presented here suggests that AEC is a good 
choice. 



394 
395 
396 
397 
398 
399 
400 
401 
402 
403 



6.3. Future work and limitations 

Despite the limitations of the set of ecosystem network models (e.g. relatively few model authors, only 
35 distinct ecosystems, some very small models n < 5, mass-balanced assumption), the consistency of the 
results suggests they are robust to model error. To verify the generality of these results, this data set 
can be exten ded to include biogeo chemically-based networks that may reflect different aspects of ecosystem 



organization (|Borrett et all 120101) as more models become available. Our preliminary analysis suggest that 



these results may hold true for biogeochemically-based networks, but this remains a testable hypothesis. 

One important direction for future research is clarifying the relationship between alternative m e thods 
for quantifying the "importance" and "effects" species have on one another ( Ulano wicz and Puccial. Il990t 
Jordan et al. . 2003 : Schramski et al. . 20061 Bauer et al.l . l201oT) . For example, Schramski et al. 's ( 2006 ') mea- 
sure of distributed control among species in a network is calculated using the input and output oriented 
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404 integral flow matrices. There arc some clear conceptual differences, but since the underlying data and an- 

405 alytical goals are the same we expect there to be an interesting relationship between AEC and distributed 

406 control. Another possible comparison would be to calculate WI beyond path lengths of 10 to see how longer 

407 indirect relationships effect this centrality metric in an ENA framework. 

408 6.4- Summary 

409 This work makes three core contributions to ecology. First, we introduce a formal method called environ 

410 centrality to quantify the relative importance of species in generating ecosystem activity. This metric builds 
4n on the centrality concept from the social sciences and Patten's environ concept, and can be considered in 

412 a manner similar to rank- abundance curves familiar to many ecologists. The main ecological advantages of 

413 environ centrality are that it is the only metric to date to include all direct and indirect effects of a weighted 

414 ecosystem network model, it is simple to calculate in the ENA framework, and it should easily apply to 

415 network models of flow in other types of complex systems. Further, it is sensitive to both extensive and 

416 intensive changes to direct flow intensities, and does not ignore the sometimes important transient dynamics 
41? like eigenvector centrality does. Second, we applied environ centrality to find evidence that support the 

418 hypothesis that ecosystems are comprised of both a few functionally dominant species and a larger group of 

419 species with more similar importance. Our results expose the central importance of detritus and decomposer 

420 species like bacteria and fungi in generating carbon throughflow activity in ecosystems. This result is not 

421 surprising, but it supports our claim that environ centrality is a useful ecological measure. Third, this 

422 work shows that indirect effects tend to homogenize the functional importance of species. This is further 

423 evidence of the important consequences that different types of indirect effects can have in complex systems 

424 like ecosystems. 
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Table 1: Formulas and references for centrality indices compared in this paper 



Metric 



Symbol Formula 



Citation 



(E^i/« + E"=i/ii)/2 

2 ^ hj k i/hj k 



Degree Centrality DC 

Weighted Degree WD 

Bctwccnncss Centrality^ BC 

Weighted Importance WI 

Eigenvector Centrality"'' EVC 

Environ Centrality AEC 



j<k;i^j (n-l)(n-2) 
Z^fc=l m 



(v + w)/2 
{EC) 



EC out 



)/2 



Freeman (1979) 
Freeman (1979) 



hjk Freeman (1979) 

Jordan et al. (2003V . Bauer et al. (2010) 
Bonacich (1972) 



t hjki represents the number of shortest paths between node j and k that pass through node i. ' » and w are the right and 



left hand eigenvectors of the dominant eigen value respectively. 
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Table 2: Fifty empirically-based trophic ecosystem network models. 



Model 


units 






nt 


ct 


T5T t 






Source 


Lake Findley 


gC m~ 2 


yr- 1 




4 


0.38 


51 


0.30 


Richcy ct al. 


(19781 


Mirror Lake 


gC m~ 2 


yr" 1 




5 


0.36 


218 


0.32 


Richev et al. 


(19781 


Lake Wingra 


gC m~ 2 


yr" 1 




5 


0.40 


1,517 


0.40 


Richev et al 


(19781 


Marion Lake 


gC m" 2 


yr" 1 




5 


0.36 


243 


0.31 


Richey et al 


(19781 


Cone Springs 


kcal m - 


2 yr" 




5 


0.32 


30,626 


0.09 


Tilly 


(19681 


Silver Springs 


kcal m - 


2 yr - 




5 


0.28 


29,175 


0.00 


Odum 


(19571 


English Channel 


kcal m - 


2 yr" 




6 


0.25 


2,280 


0.00 


J3r^linskv 


(19721 


Oyster Reef 


kcal m - 


9 

yr 




6 


0.33 


84 


0.11 


Dame and Patten 


(19811 


Somme Estuary 


mgC m~ 


' 2 d~ 




9 


0.30 


2,035 


0.14 


Rybarczyk et al 


(2003) 


Bothnian Bay 


gC m -2 


yr" 1 




12 


0.22 


130 


0.18 


Sandbcrg ct al 


(2000) 


Bothnian Sea 


gC m -2 


yr" 1 




12 


0.24 


458 


0.27 


Sandbcrg ct al 


(2000) 


Ythan Estuary 


gC m~ 2 


yr" 1 




13 


0.23 


4,181 


0.24 


Baird and Milne (1981) 


Baltic Sea 


mgC m~ 


' 2 d~ 


1 


15 


0.17 


1,974 


0.13 


Baird ct al 


(1991) 


Ems Estuary 


mgC m~ 


' 2 d~ 




15 


0.19 


1,019 


0.32 


Baird ct al 


(1991) 


Swarkops Estuary 


mgC m~ 


' 2 d~ 




15 


0.17 


13,996 


0.47 


Baird et al. (19911 


Southern Benguela Upwclling 


mgC m~ 


' 2 d" 




16 


0.23 


1,774 


0.19 


Baird et al. (19911 


Peruvian Upwelling 


mgC nT 


' 2 d~ 




16 


0.22 


33,496 


0.04 


Baird et al. 


(1991) 


Crystal River (control) 


mgC m" 


' 2 d~ 




21 


0.19 


15,063 


0.07 


Ulanowicz 


(19861 


Crystal River (thermal) 


mgC m~ 


' 2 d~ 




21 


0.14 


12,032 


0.09 


Ulanowicz 


(19861 


Charca de Maspalomas Lagoon 


mgC m~ 


~ 2 d" 




21 


0.13 


6,010,331 


0.18 


Almunia et al 


(19991 
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Baird et al. (2004b1 


Gulf of Maine 


g ww m 
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0.35 
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0.15 


Link et al 


(2008) 


Georges Bank 


g ww m 
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Link et al 


(2008) 
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32 
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Link ct al 


(2008) 
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0.51 


Monaco and Ulanowicz 


(19971 
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_ 
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Link ct al 


(2008) 
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St. Marks Seagrass, site 2 (Jan) 
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0.09 
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(1998) 
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59 
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1,353,406 


0.09 


Baird ct al. 


'2004a) 


Graminoids (wet) 


gC m~ 2 
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66 


0.18 


13,677 


0.02 


Ulanowicz ct al 


(2000) 


Graminoids (dry) 
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66 
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0.04 


Ulanowicz ct al. 


(2000) 
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68 
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0.04 


Ulanowicz et al. 
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0.04 


Ulanowicz ct al. 


(19971 


Lake Oneida (pre-ZM) 
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74 


0.22 


1,638 


< 0.01 


Miehls et al. 
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Lake Quinte (pre-ZM) 
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< 0.01 


Miehls et al. 
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< 0.01 


Miehls et al. 
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0.21 
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0.01 


Miehls et al. 


2009b1 
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94 


0.15 


3,272 


0.10 


Ulanowicz ct al. 


(19991 


Mangroves (dry) 
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yr- 1 




94 


0.15 


3,266 


0.10 


Ulanowicz ct al. ( 1999) 


Florida Bay (wet) 


mgC m~ 


" 2 yr - 


-i 
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0.12 


2,721 


0.14 


Ulanowicz ct al. 


(1998) 


Florida Bay (dry) 


mgC m~ 


-2 yr - 


-i 
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0.13 


1,779 


0.08 


Ulanowicz ct al. 


(1998) 



' n is the number of nodes in the network mode 1, C = L/n 2 is the model conncctance when L is the number of direct links or 
energy— matter transfers, TST — ^2^2 fij z i 1S the total system throughflow, and FCI is the Finn Cycling Index I IFinnl . [19801) . 
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Tabic 3: Dominate species and functional variation in 50 ecosystem network models 



Model 


CV(AEC) 


CV(R-AEC) T 


# Dominants 


Dominant Identity 


Lake Findlcy 


0.48 


0.48 







Mirror Lake 


0.51 


0.17 


2 


Sediment, Benthos 


Lake Wingra 


0.43 


0.43 







Marion Lake 


0.53 


0.14 


2 


Sediment, Benthos 


Cone Springs 


0.23 


0.23 







Silver Springs 


0.11 


0.11 







English Channel 


0.16 


0.16 







Oyster Reef 


0.33 


0.33 







Sommc Estuary 


0.58 


0.18 


1 


POM 


Bothnian Bay 


0.34 


0.34 







Bothnian Sea 


0.46 


0.46 







Ythan Estuary 


0.50 


0.50 







Baltic Sea 


0.40 


0.40 







Ems Estuary 


1.03 


0.31 


1 


Sediment POC 


Swarkops Estuary 


0.97 


0.30 


1 


Sediment POC 


Southern Bcngucla Upwclling 


0.58 


0.26 


1 


Suspended POC 


Peruvian Upwclling 


0.41 


0.41 







Crystal River (control) 


0.56 


0.27 


1 


Detritus 


Crystal River (thermal) 


0.61 


0.30 


1 


Detritus 


Charca dc Maspalomas Lagoon 


0.70 


0.40 


1 


Sedimented POC 


Northern Bcngucla Upwclling 


0.58 


0.23 


1 


POC 


Neuse Estuary (early summer 1997) 


0.79 


0.43 


2 


Sediment POC, Sediment bacteria 


Neuse Estuary (late summer 1997) 


0.83 


0.37 


2 


Sediment POC, Sediment bacteria 


Neusc Estuary (early summer 1998) 


0.80 


0.44 


2 


Sediment POC, Sediment bacteria 


Neuse Estuary (late summer 1998) 


0.81 


0.41 


2 


Sediment POC, Sediment bacteria 


Pnlf nf A/To i no 


65 


41 


2 


Ul, LI I L LIS / i KJ , KJ LI1CI IVldd UUCIlLIIUci 


Georges Bank 


68 


41 


2 


JJC LI I L LIS / i V_y v_- , KJ LIICI IVIciCI U UCIILIIUS 


Middle Atlantic Bight 


0.73 


0.45 


2 


Dctritus/POC, Other Macrobenthos 


!N ar r agansct t t3ay 


1.34 


0.39 


3 


Dntritnc; /POP Snrlimrnt POP Rartnria 

Lj U L I 1 L LlO / IT K^J v_' , O LU I11H_>II L IT KJ\^r UaLLLI Id, 










Pelagic Bacteria 


Southern New England Bight 


0.68 


0.49 


1 


Detritus/POC 


Chesapeake Bay 


1.10 


0.36 


2 


Sedimented POC, Bacteria in Sediment 










(POC) 


St. Marks Scagrass, site 1 (Jan.) 


1.00 


0.38 


2 


Sediment POC, Bcnthic Bacteria 


St. Marks Seagrass, site 1 (Feb.) 


0.94 


0.42 


2 


Sediment POC, Bcnthic Bacteria 


St. Marks Scagrass, site 2 (Jan.) 


0.86 


0.34 


2 


Sediment POC, Bcnthic Bacteria 


St. Marks Seagrass, site 2 (Feb.) 


0.86 


0.35 


2 


Sediment POC, Bcnthic Bacteria 


St. Marks Seagrass, site 3 (Jan.) 


0.60 


0.27 


1 


Sediment POC 


St. Marks Seagrass, site 4 (Feb.) 


0.75 


0.34 


1 


Sediment POC 


Sylt-R0m0 Bight 


1.17 


0.44 


1 


Sediment POC 


Graminoids (wet) 


0.98 


0.18 


2 


Sediment Carbon, Refractory Detritus 


Graminoids (dry) 


1.00 


0.28 


2 


Sediment Carbon, Refractory Detritus 


Cypress (wet) 


0.89 


0.49 


3 


Liable Detritus, Living Sediment, Verte- 










brate Detritus 


Cypress (dry) 


0.85 


0.50 


3 


Vertebrate Detritus, Liable Detritus, Liv- 








ing Sediment 


Lake Oneida (prc-ZM) 


0.49 


0.49 







Lake Quinte (pre-ZM) 


0.62 


0.28 


1 


Sedimented Detritus 


Lake Oneida (post-ZM) 


0.49 


0.49 







Lake Quinte (post-ZM) 


0.78 


0.29 


1 


Sedimented Detritus 


Mangroves (wet) 


1.17 


0.35 


3 


Sediment Carbon, Sediment bacteria, POC 


Mangroves (dry) 


1.19 


0.35 


3 


Sediment Carbon, Sediment bacteria, POC 


Florida Bay (wet) 


1.43 


0.46 


3 


Water POC, Bcnthic POC, Water Flagel- 
lates 


Florida Bay (dry) 


1.40 


0.36 


3 


Water POC, Bcnthic POC, Water Flagcl- 



latcs 



t CV(R-AEC) represents the CV of AEC scores for all species remaining after dominates have been removed. 
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Figure 1: Rank AEC curve of the Georges Bank Ecosystem network model. The grey line indicates the potential even 
distribution. Notice that the area under the curve sums to one. 
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Figure 2: Example of average environ centralities (AEC) in four realizations of a hypothetical ecosystem model. (A) Model 
structure, (B) four mass balanced realizations of the direct flow intensity matrices, G, for the model, (C) corresponding integral 
flow matrices, N, (D) AEC distributions, and (E) average eigenvector centralities AEVC. Notice that flow is from column j to 
row i in the matrices. When column sums of G do not equal 1, the remaining portion of flow is lost as exports or respiration. 
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Figure 4: Comparison of the average environ centrality (AEC) and average eigenvector centralities (AEVC) for the (A) Oyster 
Reef (B) Bothnian Bay (C) Chesapeake Bay and (D) Northern Benguela ecosystem models. The nodes have been ordered in 
ascending order of the average environ centrality to facilitate comparison. The AEC distributions were multiplied by -1 for 
plotting convenience. 
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Figure 5: Comparing the direct and integral (direct, boundary, and indirect) AEC distributions. The gray line represents the 
value if each node had the same centrality. 
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Figure 6: CV ratios quantifying the homogenization of roles in AEC when compared to AEC ect . Bars surpassing the 
vertical line at 1 represent ecosystems which exhibit centrality homogenization from indirect effects. 
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